Nonlinear collisional absorption of laser light in dense, strongly coupled plasmas 
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We present a new theoretical approach for collisional absorption of laser energy in dense plasmas 
which accommodates arbitrary frequencies and high intensities of the laser field. We establish 
a connection between laser absorption by inverse Bremstrahlung and the stopping power. This 
relation is then applied to include strong correlations beyond the mean field approach. The results 
show an excellent agreement with molecular dynamics simulations up to very high coupling strength. 
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Understanding the interaction of intense radiation with 
strongly coupled plasmas is crucial for the design and 
critical evaluation of targets for inertial confinement fu- 
sion. To accommodate the symmetry conditions, the ab- 
sorption of laser energy must be carefully determined 
starting from the early stages [2, Similar data are 
required for fast ignition by ultra-intense lasers due 
to creation of a dense plasma by the nanosecond pre- 
pulse 3|. Least understood are laser-plasma interac- 
tions that involve nonideal and partially degenerate elec- 
trons with T = (e 2 /fc B r e )(47rn e /3) 1/3 ~ 1 and n e Kl = 
n e (2irh 2 /m e kBT e ) 3 / 2 ~ 1, respectively. Such conditions 
are also created in contemporary warm dense matter ex- 
periments 0, Q and laser-cluster interactions @, 0] • 

The dominant absorption mechanism for lasers with 
intermediate intensities typical for inertial confinement 
fusion is the inverse bremsstrahlung. A first description 
of this process was presented by Dawson & Oberman 
for weak fields |8| which was later extended to arbitrary 
field strengths [9( . Due to the classical description, these 
results become inapplicable for dense, strongly coupled 
plasmas. This break-down is avoided by a rigorous quan- 
tum kinetic description applying the Green's function for- 
malism T3, U | or the quantum Vlasov equation [l2| . For 
weak laser fields that allow for a formulation within lin- 
ear response, strong electron-ion collisions were included 
into a quantum description [13j, [14 1 . 

All approaches mentioned above are formulated in 
the high-frequency limit which requires the number of 
electron-ion collisions per laser cycle to be small. In this 
limit, the electron-ion interaction has a collective rather 
than a binary character and the laser energy is coupled 
into the plasmas via the induced polarization current. On 
the other hand, binary collisions dominate laser absorp- 
tion in the low-frequency limit and a Drude-like formu- 
lation follows. At intermediate frequencies, both strong 
binary collisions and collective phenomena have to be 
considered simultaneously. Interestingly, such conditions 
occur for moderate heating at the critical density of com- 
mon Ny:Yag lasers. The strong restrictions on the ap- 
plicability of the theories above are also related to the 
assumption of nearly equilibrium electron distributions 



either in the reference frame of the ions for low frequen- 
cies or in the frame of freely moving electrons for high 
frequencies of the laser. 

In this Letter, we present a description of collisional 
absorption that bridges between the high- and low- 
frequency limits and incorporates weak collective inter- 
actions and strong binary collisions simultaneously. To 
this end, we formally split the interactions into a weak 
interaction part and additional contributions due to hard 
collisions. The latter constitutes a friction between the 
electron and ion fluids and can be treated as the stop- 
ping power of ions in the electron fluid. Thus, two basic 
energy absorption mechanism, the stopping power and 
collisional absorption, are connected which allows one to 
apply well developed models for the stopping power (see, 
e.g., Refs. HI, lij]) to the problem of laser absorption in 
plasmas. Restrictions with respect to the laser frequency 
are avoided by expanding the electrons distribution in 
a generalized Kramers-Hennenberger frame that follows 
the center of the electron fluid and is determined by the 
driving field and the friction between the two species. 
The description of the collective electron response can 
be kept almost unchanged from earlier approaches [l2| . 
Due to the use of quantum approaches for the stopping 
power and a quantum dielectric function, no ad hoc cut- 
offs must be introduced and the theory stays reliable for 
strong electron-ion interactions and degenerate electrons. 
The few assumptions made are justified by the unprece- 
dented agreement with molecular dynamic (MD) simula- 
tions [3, |2(| up to very high coupling strengths. 

Collisional absorption of laser energy is commonly 
characterized in terms of a frequency-dependent electron- 
ion collision frequency defined as [2l| 
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where the brackets denote the average over one period of 
the laser field with frequency loq and lu p = (e 2 n e /m e ) 1 / 2 
is the plasma frequency of the target electrons. The col- 
lision frequency v e i contains the same information as the 
dynamic conductivity a and is determined by the elec- 
tron current j. This current, in turn, is the first moment 
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of the electron distribution f e (p,t). We, therefore, base 
our statistical description on a general quantum kinetic 
equation for the electrons in an external field that has, 
for the homogenous plasmas considered, the form [l(| 



Ofe(p,t) 
dt 



- eE ext (t) ■ dfe{ ^ t] = I ee (p, t) + I ei (p, t) . (2) 
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Here, we will assume harmonic external electric fields: 
E ex t(t) = Eo sin (u>at). The form of the collision integrals 
on the r.h.s. is not specified at this moment. 

The current is obtained by multiplying Eq. (T5J) by 
ep/m and integrating over the free momentum. The 
terms on the l.h.s. are easily evaluated. The electron- 
electron collision term vanishes. The electron-ion colli- 
sion integral I e i is formally split into a part given by the 
dynamically screened first Born approximation and the 
remainder that contains the corrections due to hard col- 
lisions. Using the fact that the treatments based on the 
quantum Vlasov equation and the dynamic Born approx- 
imation are equivalent (compare Refs. [Hi [13]), the weak 
coupling part of the collision integral can be written as a 
self-consistent polarization P associated with the Hartree 
term. We then obtain for the electron current 
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where the collision integral I^f contains the corrections 
due to hard collisions only. 

Let us now consider the last term associated with the 
hard collisions in more detail. Since the electron quiver 
velocity is practically always much larger than the ion 
thermal motion, the ions can be considered to move with 
a common velocity V with respect to the electrons. In 
such a situation, the form of the last term in Eq. ([3]) 
is by definition the stopping power, i.e., the energy loss 
per unit length, for the ions in an electron gas (see, e.g., 
Ref. [151]). Spatial correlations in the ion component can 
be modelled as a correlated beam [171 Il8| ] . Since the last 
term contains only the contribution of hard collisions, 
it can be calculated as the difference of the full stop- 
ping power (including hard collisions) and the stopping 
power calculated from the (quantum) Lenard-Balescu 
equation describing the weak interactions. The dynam- 
ics of screening is of minor importance when describing 
the hard collisions, particularly in the low velocity case, 
but cannot always be neglected. For that reason, only 
the contribution of hard collisions is treated in this aver- 
age manner and the full dynamics is kept in the leading, 
weak coupling term of the Born series. 

It is now convenient to cast the effect of the stopping 
power associated with the hard collisions as a friction 
force or a collision frequency 
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where x and V are pointing in the same direction as the 
external field Eo and V is the magnitude of the time- 
dependent velocity between the electron and ion fluids. 
Again, the stopping power and the generalized friction 
coefficient R contain the effect of hard collisions only. 

In general, the above model must be solved numeri- 
cally. However, if the friction coefficient R is not velocity- 
dependent as during the linear increase of the stopping 
power at small V, one can proceed analytically. This 
regime is determined by the opposing forces created by 
the external field and friction due to hard collisions. As 
the stopping power is linear in V up to roughly the ther- 
mal velocity of the electrons, we have 
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with the thermal velocity v t h = (^B^e/me) 1 / 2 and the 
free quiver velocity vq = eEo/muJo. In the high frequency 
limit, we recover v < v t h as the restriction for the field 
strength while vq < [1 + (v^/bj )] v th follows in the low- 
frequency, i.e. Vgf/uiQ 3> 1, limit. Thus, the region where 
R = const can be applied is extended to significantly 
higher field amplitudes for strongly coupled plasmas and 
low to intermediate frequencies. 

Given that V < V max , we can formally solve Eq. ([3]) 
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where we introduced 7 = [1 + (^ c /w ) 2 ] -1 . The total 
current consists of contributions from the moving carriers 
and the polarization. The strong collisions result in the 
second term in the brackets which is in phase with the 
driver field. If the contribution from the polarization 
can be neglected as in the low- frequency limit, Eq. ([3]) 
gives the well-known Drude formula for the conductivity: 
ctd(wo) = Wp/4tt(iA c - iuj ). 

By inserting Eq. ([6|) into the definition of the collision 
frequency (fT]), one obtains 



v el (u) = 71/, 



(7) 

The contribution of the hard collisions in the first term 
has the form of the real part of the Drude conductivity. 

Deriving an explicit expression for the second term 
in Eq. ([7]) associated with the polarization is mainly 
straightforward. Differences to earlier works arise from 
the additional hard collisions which define a generalized 
frame for the center of mass of the electron fluid. Here, 
we only sketch the derivation using classical arguments; 
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the quan tum version of the derivation follows the lines of 
Ref. [l2| and will be published elsewhere [22j. 

To find the polarization field, one uses its relation to 
the self-consistent potential of the plasma: P = — V<I>. 
The field is, in turn, determined by the Poisson equation 



V 2 <& = 4nn e e J dp f e (p) - AirZe'^8{r - n) , 



(8) 



where we consider point-like ions positioned at tv A 
direct solution like for the Vlasov equation is not possible 
here as the collision integral depends on the momentum 
p. If we, however, keep in mind that the hard collisions 
are treated as an average friction force acting on all ions 
and electrons, we get a Vlasov-like equation of the form 



dfefat) 
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which is consistent with the current balance equation ([3J . 
Such a treatment implies that all individual electrons feel 
the same, average friction force as the electron fluid. 

Next a frame where the electron distribution can be 
expanded about its equilibrium form has to be defined. 
Due to the incorporation of hard collisions, it differs from 
the one of freely quivering electrons and is determined 
by the external field and the friction force (generalized 
Kramers- Hennenberger frame). The transformation to 
this frame is given by 
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where the abbreviation e = — eEo/mtjQ is used. In this 
frame, the electron distribution can be linearized about 
its equilibrium /o which yields 



dfi(m e u,t) e d$i(p,t) df (m e u,t) 



dt 
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The equation for the lineraised potential <E>i is easily ob- 
tained from the Poisson equation ||5J) using the first order 
electron fluctuation fi and a coordinate translation ac- 
cording to (|10p . After Fourier transformation, the poten- 
tial can be written in the form $i(fe,o;) = S(fe, uS)/e(k, u) 
where e(k,u) is the dielectric function. The background 
source term £(fe,u;) can be expanded using Bessel func- 
tions of the first kind to yield 

xS(uj + [n + m](j ) ^ ^ ' ( 12 ) 

j 

where the arguments of the Bessel-functions are given by 
£ = 7&- e and £ = 7iA c fc- e/u>Q. The time-dependent fields 



follow by inverse Fourier transformation and one obtains 
P(k,t) = ik&i(k, t) for the polarization field. This field is 
now transformed back to the rest frame of the ions and 
averaged over ion positions. In the last step, the sum 
over exponentials in Eq. fT3]) becomes the static ion-ion 
structure factor Su(k). For the polarization in coordinate 
space, one gets 

P{r,t) = i J dk k cj>(k,t) Sn(k) 

x exp(i£\sm(u t) + — cos(w t)l") . (13) 



This field can be easily averaged over a laser cycle and 
then inserted into the expression ((7| for the collision fre- 
quency which yields 
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where a — (iVQu; p /2TT 2 )(Ze 2 /m e v^). Clearly, the known 
limiting cases can be readily retrieved. The last term 
in line 3 dominates for small laser frequencies giving a 
Drude-like expression. The results from quantum Vlasov 
and Born approximation follow in the weak coupling limit 
with fgf— ->0. The result of Decker et al. 9] follows also 
from expression (fT4"|) for — > if the classical dielectric 
function is used. In this case, the integral must be trun- 
cated at k m ax to avoid the divergence at small impact 
parameters. 

The form of expression (|14[) should also clarify why the 
collision term was split into a weak coupling and a hard 
collision part: only a Drude-like term would have been 
obtained if the total electron-ion collision integral was 
evaluated by using its relation to the stopping power. 
Moreover, nonlinear contributions associated with the 
higher order Bessel functions would be neglected. Here, 
the full dynamics in the leading order term is kept and 
the effects of strong collisions are incorporated as well. 

The derivation above was sketched classically, but can 
be done quantum-mechanically as well. The only dif- 
ference that arises is the form of the dielectric function 
e(k, uj) which will be classical or quantum RPA (compare 
Ref. Q with Refs. [HI, [13]). To obtain the results shown 
below in Fig. [T] we used the quantum form including 
degeneracy corrections for e(k,uj). The stopping power, 
that defines the frequency of hard collisions , is also 
obtained from a quantum kinetic descriptions. As a re- 
sults all integrals can be performed to infinity and no ad 
hoc cutoffs must be introduced. 

The main difference of expression (|14p to other results 
for the collision frequency is the incorporation of hard 
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Polarization (14) 
Hard collisions (14)' 
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of other cutoffs). For a coupling strength of T w 1, the 
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quantum theories of Refs. [ll|, [l2| also start to disagree 
with the simulation data. This behavior can be traced 
back to i) the neglect of hard collisions and ii) the use of 
a freely moving reference frame. As both shortcomings 
are overcome in our approach we find excellent agree- 
ment with data from MD simulations up to high cou- 
pling strength. Plotting both contributions of Eq. ([M]) 
separately reveals that the Drude term (hard collisions) 
dominates for high coupling strengths and defines here 
the shape of the curve. This adds another nonlinearity 
to the result for more highly charged ions since the rele- 
vant stopping power is here not proportional to Z 2 flij ]. 

In conclusion, a new approach for collisional absorption 
of laser light in dense plasmas was presented. It incorpo- 
rates the dynamic response and hard electron-ion colli- 
sions. The latter were introduced by a newly established 
connections to the stopping power which also allowed one 
to define a reference frame that includes both field effects 
and friction. Although only results for the quasi-linear 
regime vo/vth < 1 were presented, this approach is not 
limited by this requirement and can be easily extended 
to higher field amplitudes, correlated ions, and multiple 
ionization stages. 

The authors thank J. Vorberger (CFSA, Warwick) for 
fruitful discussions and EPSRC for financial support. 
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FIG. 1: Collision frequency v e i versus coupling parameter T 
for a hydrogen plasma with n e = 10 22 cm -3 and a laser field 
with u>o/ui p = 3 and vo/v t h =0.2. Solid line: Eq. (|14[1 ; punc- 
tured lines: contributions of hard collisions and polarization 
to Eq. (|14|l : dashed line: quantum results as in Refs.Jll"l.ll2T|; 
dash-dotted line: classical results of Decker et al. 9]) with 
an integral cut off at k max =m e w 2 h /Ze 2 ; symbols: results of 
numerical simulations. 



electron-ion collisions. This has been done by introduc- 
ing a general friction forces related to the stopping power 
of the ions in an electron gas. Many models for the stop- 
ping power have been developed [25| . few include hard 
collisions. Within quantum statistical theory, they can 
be described by a T-matrix approach based on the quan- 
tum Boltzmann equation. The related cross sections are 
then calculated by numerical solutions of the Schrodinger 
equation (IBj . Dynamic screening effects can be added 
applying the Gould-DeWitt scheme 24J or by velocity- 
dependent screening length [lq |. While both methods 
coincide for weak and intermediate coupling, only the 
latter agrees with simulation data for very high beam- 
plasma coupling [3, [2^| . The low velocity part of this 
model can be expressed in an analytic fit of the form 



logHth.u.]) = exp(-0.0735x 2 +1.3373x- 1.8511) , (15) 

where x = ZT 3 ' 2 and the total collision frequency v is 
given in thermal units. To obtain the contribution of 
hard electron- ion collisions needed in expression (fT4|) . the 
simple Lenard-Balescu form must be subtracted. 

In Fig. [1] the results of our new approach are com- 
pared to other theories and simulation data 1^, IS 23j |. 
According to the conditions considered, the stopping 
power was here calculated using the combined model of 



Refs. [15|, [24j, i.e., is given by the difference of the 
T-matrix and the static Born terms. The ions can be 
treated individually since Su ~ 1 even for the highest 
coupling strength considered. As expected, all theories 
agree for weakly coupled plasmas, but large deviation 
occur for strong coupling. The classical description is 
clearly not applicable here (see Ref. [22| for a discussion 
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